Supersymmetric Electroweak Phase Transition: 
Baryogenesis versus Experimental Constraints 



James M. Cline, Guy D. Moore 

Physics Department, McGill University, 3600 University Street, Montreal, Quebec, Canada H3A 2T8. 

We use the two loop effective potential to study the first order electroweak phase transition in 
the minimial supersymmetric standard model. A global search of the parameter space is made to 
identify parameters compatible with electroweak baryogenesis. We improve on previous such studies 
by fully incorporating squark and Higgs boson mixing, by using the latest experimental constraints, 
and by computing the latent heat and the sphaleron rate. We find the constraints tan/3 > 2.1, 
TTij < f 72 GeV, and ruh < 107 GeV (becoming more or less restrictive if the heavy stop has mass 
less than or greater than 1 TeV). We also find that the change in tan/3 in the bubble wall is rarely 
greater than 10"'', which severely constrains the mechanism of baryogenesis. 
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The niost predictive theory of the origin of baryonic 
matter at present is electroweak baryogenesis, since it 
relies on assumptions about new physics that is being 
searched for at LEP and Fermilab. For sufficient baryo- 
genesis, the electroweak phase transition (EWPT) must 
be strong; the requirement is roughly v/Tc > 1, with v 
the Higgs vacuum expectation value and Tc the critical 
temperature. This is unachievable in the standard model. 
One of the main contenders for the new physics needed is 
supersymmetry, whose simplest manifestation is the Min- 
imal Supersymmetric Standard Model (MSSM). Much 
progress has been made in the last few years toward iden- 
tifying which regions of the MSSM parameter space are 
compatible with electroweak baryogenesis It is 

agreed that a light Higgs boson and a light top squark are 
needed, but estimates vary as to just how light. One goal 
of the present work is to refine these estimates within the 
framework of the effective potential approach. 

Computing the properties of the EWPT is subtle be- 
cause the power-counting arguments for perturbation 
theory at zero temperature are modified at high temper- 
ature, such that the loop expansion for the free energy 
(effective potential) must be resummed. Even after re- 
summing, the effective potential (EP) is not guaranteed 
to be reliable at small values of v, due to infrared diver- 
gences coming from the small masses of the transverse 
gauge bosons and the light Higgs boson. There are three 
principal ways of dealing with this problem: (1) com- 
pute phase transition properties fully numerically, on the 
lattice (2) integrate out the heavy, nonlR-divergent 
modes analytically ("dimensional reduction") and save 
the lattice for the effective theory of the dangerous light 
modes [^-0]; and (3) improve the effective potential by 
computing to higher order in the loop expansion |^,^,|l2| . 

The first approach, because of its high numerical cost, 
is impractical from the standpoint of exploring large re- 
gions of parameter space. Dimensional reduction is bet- 
ter in this respect because it maps the full MSSM param- 
eter space onto a much smaller set of parameters in the 
effective high-T theory for the light modes, thereby re- 



ducing the number of lattice simulations needed. But DR 
has its limitations, especially when it comes to the effects 
of superrenormalizable interactions which can induce a 
large number of unsuppressed high-dimension operators 
in the effective theory to be studied on the lattice. The 
At and fi parameters of the MSSM come with precisely 
these kinds of interactions. 

On the face of it, computing the effective potential 
to higher order in perturbation theory would not seem 
promising since the convergence of the perturbation se- 
ries is supposed to be poor. However experience with the 
standard model shows that the two-loop EP works quite 
well. In the MSSM there are also indications that the 
two-loop EP gives converging results, since they are in 
fairly good agreement with recent lattice computations: 
the lattice gives values of v /T (a measure of the strength 
of the transition) which are about 10% higher than those 
of the EP. This, together with the ability it affords for 
quickly combing the full MSSM parameter space, makes 
it worthwhile to investigate the two-loop EP. Further- 
more, detailed properties of the phase transition like the 
nucleation temperature Tnuc, sphaleron energy inside the 
bubbles, and the variation of the two Higgs field VEV's 
inside the bubble walls, are much more readily available 
from the EP than from lattice studies. 

The important contributions to the two-loop EP in the 
MSSM have been calculated by references ||,|,|l|. The 
first reference included squark mixing, while assuming 
the heavy Higgs bosons are decoupled. The third allowed 
for lighter Higgs bosons, but ignored squark mixing. We 
have generalized these results to incorporate both effects. 
This is desirable because most baryogenesis mechanisms 
rely upon tan/3, the ratio {v2/vi) of the two Higgs VEV's, 
changing from the interior to the exterior of the bubbles, 
which is only possible if the heavy Higgs bosons are not 
decoupled. Squark mixing is also expected because the 
phase (and magnitude) of the fi parameter must not be 
too small, since this is the principal source of CP violation 
p^- Q. The fj, parameter appears in the off-diagonal 
term of the stop mass matrix, while the diagonal term 
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for the right-handed stop must be small to get a strong 
enough phase transition. Thus significant squark mixing 
is a possibility, which in fact is realized (figure 2 below). 

A recent development in the literature is to consider 
the formation of squark condensates at the beginning of 
a two-stage transition, where color and electric charge 
are temporarily broken [^,0. To recover our known 
world, the color broken state must later copiously nucle- 
ate bubbles of the conventional SU(2)-breaking ground 
state of the standard model. If this could happen, it 
would considerably strengthen the EWPT. However, this 
nucleation process is heavily suppressed iQ, leading to 
the same problem that killed "old" inflation-the CCB 
vacua would inflate and vastly dominate the spatial vol- 
ume of the universe. Hence, at all temperatures above 
the nucleation temperature for the electroweak broken 
phase, the effective potential for the right stop must have 
a stable or sufficiently metastable symmetric minimum 
to prevent stop condensation. Analyzing the one loop 
effective potential for the right stop, we find that nu- 
cleation of the color broken vacuum happens whenever 
(T) < 0.035 g^T^. We therefore discard parameters 
where this happens at any T > T^uc- 

The most relevant laboratory constraints concern mh, 
mi and (the lightest Higgs boson, light stop and sbot- 
tom masses, respectively), Ap (the contributions of the 
stops and sbottoms to the p parameter), and the ex- 
clusion of charge- and color-breaking minima. The lat- 
est experimental limit on mh depends on sin^(a — /3), 
where a describes the composition of the light Higgs 
field h through h = sinaff? + cos aH^ The 95% 

c.l. excluded region is roughly given by the intersec- 
tion of rrih > (69 + 19sin^(Q; - /3)) GeV/c^ and nih > 
(76 - 11.5sin2(a - /?)) GeV/c^ This allowed region 
and our accepted Monte Carlo points are shown in figure 
1. The limit sin^(Q; — /3) — > 1 corresponds to a heavy A° 
boson and a standard-model-like Higgs sector, with only 
one light Higgs boson. The phase transition is typically 
strongest in this regime. 

For the squark masses we use the preliminary ALEPH 
limit of mj > 82 GeV, which is left-right mixing inde- 
pendent, and nij^ > 79 GeV Concerning the devi- 
ation in the p parameter; the standard model value of 
Ap (also known as ej) is already 1.5cr larger than the 
experimental value [g^ for mh ^ 100 GeV. We con- 
strain the squark contribution to Ap to be less than ap- 
proximately one additional standard deviation, namely 
Ap < 1.5 X lO^-^. Chargino/neutralino searches also con- 
strain IpI > 100 GeV |2|. 

We search this allowed parameter space for those val- 
ues that give a strong enough phase transition. Our cri- 
terion is that the integrated rate of sphaleron transitions 
since the phase transition reduces the baryon asymme- 
try by just one e-folding. Writing the sphaleron rate per 
unit volume as Ts and the sphaleron energy as Egph, the 
bound is Ts < TcHt pCtl, where 
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all evaluated at the reheat temperature Tr ^ Tc- We 
find the sphaleron energy by solving for the sphaleron 
configuration using the two loop effective potential, and 
we multiply both v/Tr and and Esph/Tr by a correction 
factor of 1.1 suggested by the lattice results 0. The 
frequencies of the MSSM input parameters which pass 
all these cuts, as well as histograms for some derived 
quantities, are shown in figure 2. 
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Figure 1: (a) Experimentally allowed values of sin^(a — /3) 
and TUh are to the right of the line; Monte Carlo generated 
points are those consistent with electroweak baryogenesis. (b) 
Relation between sin'^(a — fJ) and tua- 

tanp _ rfiu niQ m^. 
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Figure 2: Frequencies of baryogenesis-allowed parameters. 
Masses are in GeV; is the stop mixing angle. 

The strength of the phase transition depends most sen- 
sitively on tan /3 and m^, which in turn determine the 
masses of the lightest Higgs boson and the top squarks. 
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Qualitatively, the dependences can be understood as fol- 
lows. In the absence of the light stop, and ignoring 
the experimental constraint on the Higgs boson mass, a 
strong phase transition requires a small value mh and 
hence small tan/3. This is because the quartic terms 
of the tree-level potential, 5^(|ifip — |772p)^, are fiat 
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along the direction of tan/3 = 1 (|i?i| = |i?2|), so the 
effective quartic couphng A is minimized for tan/3 ^ 1, 
which helps the strength of the phase transition, since 
v/T ~ <7'^/A. The experimental bound on m/i translates 
into a lower limit on tan /?, which excludes the whole re- 
gion where baryogenesis is viable. However, this can be 
counteracted if the (mostly right-handed) stop is suffi- 
ciently light. Thus the contours of constant v/T resem- 
ble hyperbolas in the tan (3-fhu plane, as shown in figure 
3. We define fhu = mfj/\mu\ so that mjj has the same 
sign as mfj. 
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Figure 3: Contours of v/T (solid) and Higgs mass (dashed) in 
the plane of tan/3 and ihu = mfj/\mu\ (in GeV), for mg — 
100 and 500 GeV, respectively, at zero squark mixing (jj, — 
At = 0). The potential has color-breaking minima in the 
black regions near rhu = —70 GeV. 

The next most important parameter is the soft- 
breaking mass term for the left-handed top squark, mg. 
It affects the strength of the phase transition almost ex- 
clusively through a radiative correction to the Higgs self- 
coupling, proportional to ln( (mg -|- )/m^). Thus, at 
fixed tan/3, raising mg increases mh and weakens the 
phase transition. (This weakening can be compensated 
by sufhciently lowering the right-stop mass.) However 
the allowed range of tan/3 consistent with the experi- 
mental limits on m^ increases with mg, as shown in fig- 
ure 4, which has the scatter plots from the Monte Carlo 
for rrih and tan/3 as functions of mg = mg/lOO GeV. 
The upper limit on as a function of mg is simply 
the maximum theoretically allowed value in the MSSM. 
The fitted functions for the upper limit on m^j and the 
corresponding lower limit on tan (3 are 



ruh < 85.9 -t- 9.2 In(mg) GeV 
tan/3 > (0.03 + 0.076 mg- 0.0031 m|)-\ 
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generalizing the findings of ref. As for the small- 

est possible values of mh, we see from figure 1 that the 
scarcity of points with m/j < 88 GeV is due to the small 
probability of getting a strong phase transition when 
ruA < 100 GeV. 

We have somewhat painstakingly reconsidered the cri- 
terion for a strong enough phase transition. We first find 
the nucleation temperature r„, which is lower than Tc, 



and the latent heat of the phase transition. From these 
we get the temperature to which the universe reheats, T^, 
on completion of the phase transition. We compute the 
sphaleron energy using the two loop effective potential at 
this temperature, and compare the resulting sphaleron 
rate to the expansion rate of the universe, also ac- 
counting for the time dependence of F^, eq. (|l|). In figure 
5a we show the correlation between the rigorous measure 
of baryon dilution, — InF^/Fcrit, and v{Tc)/Tc. Here, in 
contrast to figure 2, Fg/Fcrit is computed without apply- 
ing the lattice correction factor to v{Tc)/Tc. The points 
below the hue — InFs/Fdit — should be discarded, ac- 
cording to eq. (|l|), and are only retained if we account 
for the fact that the effective potential underestimates 
v/T. The correlation between the correct criterion and 
v(Tc)/Tc is good but not perfect. The smallest allowed 
value of v{Tc)/Tc is 1.05, and the largest rejected value 
is about 1.15. 
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Figure 4: Correlation of ruh and tan/3 with the left-handed 
top squark mass parameter mg = mg/lOO GeV. Heavy lines 
show the approximate limiting values, eq. {W. 
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Figure 5: (a) Correlation of the sphaleron rate with v/T; 
below the line would be ruled out baryon preservation, eq. (|l|); 
(b) AT = T„- Tc in GeV versus A = (latent heat)/r"*, fit by 
the function AT = -0.5 + 0.5A - 2.9A^ 

Another measure of the strength of a first order phase 
transition is the latent heat L, which is the difference 
between phases of dV/dhiT at the critical temperature, 
where V is the effective potential. It correlates strongly 
with the amount of supercooling before bubble nucle- 
ation. Figure 5b shows the relation between T„ — Tc 
and the latent heat. The latent heat also determines the 
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reheat temperature: L is the heat available for increasing 
the plasma energy density, 
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Therefore, writing Apc 
that 



g,7r2(T4 - r,f)/30, it follows 



(T, - Tr)/{T, - T„) ~ 1 - L/Ap, 
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and the universe reheats back to if the ratio L/ /^pc^n 
ever exceeds unity. However we find that L/lS.pc^n has 
an average value of 0.29 and never falls outside the range 
[0.17, 0.42], so reheating to Tc does not ever occur. 

We have also investigated the value of tan (3 = V2 /vi = 
{H2) I (Hi) inside the bubble wall. This quantity is of in- 
terest because most (but not all electroweak baryo- 
genesis mechanisms in the MSSM predict the baryon 
asymmetry is proportional to an average of ?;29ui — i'i9w2 
over the the bubble wall, where v goes from to the 
broken phase value Vc- We characterize the variation 
of tan /3 with v by minimizing the effective potential 
in the angular direction for each fixed value of v ~ 
((ifi)^ + (i?2)^)"^^^, and computing A/3 = ma.Xv[v{f3{v) — 
l3{vc))]/vc- This definition of A/3, which differs from 
ref. |13|'s, corresponds more closely to the quantity which 
enters into computations of the baryon asymmetry. Like 
ref. we find that this quantity never exceeds 0.02, 
and is typically 10 times smaller. It is also never large if 
the A'^ boson is heavy. 
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Figure 6: (a) Maximum deviation in weighted Higgs VEV 
orientation, A/3 = ma,Xv[v{/3{v) — /3{vc))]/vc, inside bubble 
wall, as a function of m^o; (b) distribution of A/3 values. 

The most interesting issue confronting electroweak 
baryogenesis in the MSSM is whether the LEP 200 run 
will be able to rule it out. The answer seems to be "al- 
most, but not completely." The final center of mass en- 
ergy Vs = 200 GeV wiU exclude nih up to 107 GeV (H). 
However if mg — 2 TeV, we find parameters with m/i as 
high as 116 GeV yet consistent with electroweak baryo- 
genesis. Such a Higgs boson might be discovered in Run 
II at Tevatron 
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studied bubble nucleation issues and A/3 prior to this 
work. 
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